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The highest superconducting temperature T c observed in any elemental metal (Li with T c ~ 
20 K at pressure P~ 40 GPa) is shown to arise from critical (formally divergent) electron-phonon 
coupling to the transverse T\ phonon branch along intersections of Kohn anomaly surfaces with the 
Fermi surface. First principles linear response calculations of the phonon spectrum and spectral 
function a 2 F(uj) reveal (harmonic) instability already at 25 GPa. Our results imply that the fee 
phase is anharmonically stabilized in the 25-38 GPa range. 

PACS numbers: 



The recent observations of superconductivity in 
fee Li up to T c = 14 K in near-hydrostatic fee-phase 
samples, p) and as high as 20 K in non- hydrostatic 
pressure cells, 0, in the pressure range 20 GPa 
<P< 40 GPa provides almost as startling a devel- 
opment as the discovery£| in 2001 of T c = 40 K in 
MgB2. Lithium at ambient conditions, after all, is a 
simple s-electron metal showing no superconductiv- 
ity above 100 /zK.[5j What can possibly transform it 
into the best elemental superconductor known, still 
in a simple, monatomic, cubic phase? There is no 
reason to suspect a magnetic (or other unconven- 
tional) pairing mechanism, but it seems equally un- 
likely that it transforms into a very strongly coupled 
electron-phonon (EP) superconductor at readily ac- 
cessible pressures. 

The strength of EP coupling in Li has attracted at- 
tention for some time. Evaluations based on empir- 
ical pseudopotentials@ early on suggested substan- 
tial coupling strength A=0.56 and hence readily ob- 
servable superconductivity (T c > 1 K); more recent 
calculations relying on the rigid muffin-tin approx- 
imation (RM TA) reached a similar conclusion]]], @ 
and led to prediction of remarkably high T c ~ 70 K 
under pressure. Q None of these studies actually cal- 
culated phonon frequencies, relying instead on esti- 
mates of a representative phonon frequency lu based 
on the Debye temperature, which is only an extrap- 
olation from the q — » phonons. Linear response 
calculations of the phonons and EP coupling 9] in 
bec Li confirmed that superconductivity would oc- 
cur in bee Li (A = 0.45), but superconductivity is 
not observed due to the transformation into the 9R 
phase with 25% weaker coupling. Experimentally, 
superconductivity only appears above 20 GPa in the 
fee phase. 

In this paper we focus on the monatomic fee phase 
that is stable in the 20-38 GPa range. After pro- 



viding additional characterization of the previously 
discussed [13 [H [II [13 evolution of the electronic 
structure under pressure, we analyze the implica- 
tions of the Fermi surface (FS) topology for proper- 
ties of Li. To study A microscopically we focus on 
the decompositionjlTj into mode coupling strengths 
X Qu , where A = (1/3JV) ^Qv =< Aq„ > is the 
Brillouin zone (BZ) and phonon branch (v) aver- 
age. We find that increase of pressure leads to very 
strong EP coupling to a specific branch in very re- 
stricted regions of momentum space determined by 
the FS topology; these features are dire ctly analo- 
gous to the focusing of coupling strength )!! . fl5l Hq 
in MgB2. Unlike in MgB2, tuning with pressure 
leads to a vanishing harmonic frequency at ~25 GPa, 
beyond which the fee phase is stabilized by anhar- 
monic interactions. 

The volume at 35 GPa is 51% of that at P=0, so 
the conduction electron density has doubled. The 
shift in character from s to p is analogous to the 
s — > d crossover in the heavier alkali metals. |is[ 
The occupied bandwidth increases by only 14%, 
much less than the free electron value 2 2 / 3 -l = 59%; 
this discrepancy is accounted for by the 55% in- 
crease in the k=0 band mass (mf,/m=1.34 at P=0 
to m(,/m=2.08 at 35 GPa). At P=0 in the fee phase 
the FSs are significantly nonspherical and just touch 
at the L points of the BZ; necks (as in Cu), where 
the p character is strongest, grow with increasing 
pressure, and the FS at 35 GPa is shown in Fig. ^ 
colored by the Fermi velocity. The topology of the 
FS plays a crucial role in the superconductivity of 
Li, as we discuss below. 

The coupling strength A is the average of mode 
coupling constants |T7j 




FIG. 1: (color online) Top figure: Fermi surface of Li 
at 35 GPa plotted in a cube region around k=0 and 
colored by the value of the Fermi velocity. Red (belly 
areas) denotes fast electrons {v^ ax = 9xlO r cm/s), blue 
(on necks) denotes the slower electrons (u™ n = 4xl0 7 
cm/s) that are concentrated around the FS necks. The 
free electron value is 1.7x 10 8 cm/s. Middle panel: Fermi 
surfaces with relative shift of 0.71(1,1,0) (i.e. near the 
point K) indicating lines of intersection. Bottom panel: 
the light areas indicate the "hot spots" (the intersection 
of the Kohn anomaly surfaces with the Fermi surface) 
that are involved in strong nesting and strong coupling 
at Q=0.71(l,l,0) (see Fig. 0. These include the necks, 
and three inequivalent lines connecting neck regions. 



with magnitude determined by the EP matrix el- 
ements M^\ + q and the nesting function £(Q) de- 
scribing the phase space for electron-hole scattering 
across the FS (E F =0), 

ttQ) = ^J26(e k )5(e k+Q )cx £ ^ ^ (2) 
N ^ J \v k x v k+Q \ 

Here the integral is over the line of intersection of 
the FS and its image displaced by Q, v k = V k e k is 
the velocity, and N(0) is the FS density of states. 
Evidently £(Q) gets large if one of the velocities gets 
small, or if the two velocities become collinear. 

Note that ^Eq£(<3) = I N (°)] 2 ; the topology of 
the FS simply determines how the fixed number of 
scattering processes is distributed in Q. For a spher- 
ical FS £(Q) oc j^9(2k F — Q); in a lattice it is 
simply a reciprocal lattice sum of such functions. 
This simple behavior (which would hold for bcc Li 
at P=0, for example) is altered dramatically in fee 
Li, as shown in Fig. Elfor P=35 GPa (the nonphysi- 
cal and meaningless t^t divergence around T should 
be ignored). There is very fine structure in 
that demands a fine k mesh in the BZ integration, 
evidence that there is strong focusing of scattering 
processes around the K point, along the T-X line 
peaking at | r-X=Xif , and also a pair of ridges (ac- 
tually, cuts through surfaces) running in each (001) 
plane in K-X^-K-Xif-K-X^-K squares. Some addi- 
tional structures are the simple discontinuities men- 
tioned above, arising from the spherical regions of 
the FS. 

Structure in £(Q) arises where the integrand in 
Eq. [^becomes singular, i.e. when the velocities at k 
and k + Q become collinear. The FS locally is either 
parabolic or hyperbolic, and the nature of the singu- 
larity is governed by the difference surface which also 
is either parabolic or hyperbolic. In the parabolic 
case (such as two spheres touching) £(Q) has a dis- 
continuity. In the hyperbolic case, however, £(Q) di- 
verges logarithmically. Such divergent points are not 
isolated, but locally define a surface of such singular- 
ities (or discontinuities, in the parabolic case). The 
ridges and steps visible in Fig. [21 are cuts through 
these singular surfaces (more details will be pub- 
lished elsewhere) ; the intensity at K arises from tran- 
sitions from one neck to (near) another neck and is 
enhanced by the low neck velocity Roth et al. have 
pointed out related effects on the susceptibility fl9| 
(which will analogously impact the real part of the 
phonon self-energy), and Rice and Halperin[2(j have 
discussed related processes for the tungsten FS. In 
the susceptibility (and hence in the phonon renor- 
malization) only FS nesting with antiparallel veloc- 
ities gives rise to Q-dependent structure. This ex- 
plains why the ridge in along the T-X line (due 
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FIG. 2: (color online) Surface plots of the nesting func- 
tion £(Q) at 35 GPa throughout three symmetry planes: 
(010) r-X-W-K-W-X-r; (001) r-K-X-r; (HO) r-K-L-K- 
X-r. The F point lies in the back corner. The dark (red) 
regions denote high intensity, the light (blue) regions de- 
note low intensity. The maxima in these planes occur 
near K and along F-X. To obtain the fine structure a 
cubic k mesh of (27i"/a)/160 was used (2xl0 6 points in 
the BZ). 



to transitions between necks and the region between 
necks) does not cause much softening (see below); 
there will however be large values of Aq„ because its 
structure depends only on collinearity. 

Divergences of £(Q), which we relate to specific 
regions of the FS shown in the bottom panel of 
Fig- H (mostly distinct from the flattened regions 
between necks discussed elsewhere^!), specify the 
Q regions of greatest instability. However, insta- 
bilities in harmonic approximation {ujq v — > 0) may 
not correspond to physical instabilities: as the fre- 
quency softens, atomic displacements increase and 
the lattice can be stabilized to even stronger cou- 
pling (higher pressure) by anharmonic interactions. 
Thus, although we obtain a harmonic instability at 
Q^K already at 25 GPa, it is entirely feasible that 
the system is anharmonically stabilized beyond this 
pressure. We infer that indeed the regime beyond 
25 GPa is an example of anharmonically stabilized 
"high T c " superconductivity. 

The phonon energies and EP matrix elements have 
been obtained from linear response theory as im- 
plemented in Savrasov's full-potential linear muffin- 
tin orbital code.j^ Phonons are calculated at 72 
incquivalent Q points (a 12x12x12 grid), with a 
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FIG. 3: (color online) Top panel: Calculated phonon 
spectrum (interpolated smoothly between calculated 
points (solid symbols) of fee Li along the T-K direction, 
at the four pressures indicated. The T\ (lowest) branch 
becomes harmonically unstable around K just above 20 
GPa. Bottom panel: calculated spectral functions o?F 
for P = 0, 10, 20, and 35 GPa. Note that, in spite of the 
(expected) increase in the maximum phonon frequency, 
the dominant growth in weight occurs in the 10-20 meV 
region. 



40x40x40 grid for the zone integration. To illus- 
trate the evolution with pressure, we use the fee lat- 
tice constants 8.00, 7.23, 6.80, and 6.41 bohr, cor- 
responding approximately to 0, 10, 20, and 35 GPa 
respectively (and we use these pressures as labels). 

The phonon spectrum along T-X behaves fairly 
normally. The longitudinal (£) branch at X hard- 
ens from 45 meV to 87 meV in the 0-35 GPa range, 
while the transverse (T) mode at X remains at 30-35 
meV. Along T-L the behavior is somewhat more in- 
teresting: again the C branch hardens as expected, 
from 40 to 84 meV, but the T branch remains low 
at 15-17 meV at the L point and acquires a notice- 
able dip near the midpoint at 35 GPa. The impor- 
tant changes occur along the (110) T-K direction as 
shown in Fig. [3] the C and T2 branches harden con- 
ventionally, but the < 110 > polarized T\ branch 
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softens dramatically around the K point, becoming 
unstable around 25 GPa. At 35 GPa this mode is 
severely unstable in a substantial volume near the K 
point (not only along the T-K line). 

We have evaluated the EP spectral function 
a 2 F(u) using our mesh of 72 Q points and the tetra- 
hedron method. Due to the fine structure in £(Q) 
and hence in Aq„, numerically accurate results can- 
not be expected, but general trends should be ev- 
ident. The resulting spectra are displayed in Fig. 
[3Jb) for each of the four pressures, showing the hard- 
ening of the highest frequency C mode with pressure 
(43 meV — > 83 meV). The most important change is 
the growth in weight centered at 25 meV (10 GPa) 
and then decreasing to 15 meV (20 GPa) beyond 
which the instability renders any interpretation at 
35 GPa questionable. The growing strength is at low 
energy; note however that this region is approach- 
ing the energy uj op t = 27rfcsT c s=s 10 meV which 
Bergmann and Rainer 22] found from calculation of 
5T c /5a 2 F(ui) to be the optimal position to concen- 
trate the spectral weight. These a 2 F spectra give 
the values of uJi og , < lu 2 > 1 / 2 , and A given in Ta- 
ble [I] The commonly chosen value n* = 0.13 in the 
Allen-Dynes equation|2^| (which describes the large 
A regime correctly) gives observable values of T c = 
0.4-5 in the 1-10 GPa range, but Li is not fee at these 
pressures. The 20 K obtained for 20 GPa is satis- 
fy ingly close to the range of observed T c , and could 
be depressed to the observed value by anharmonic 
interactions or by a larger value of /x*. 

We have shown how Fermi surface topology can 
concentrate scattering processes into specific sur- 
faces in Q-space, and even in alkali metals can lead 
to very strong coupling to phonons with these mo- 
menta, and can readily drive lattice instability. To 
enhance A, it is necessary in addition that the large 
regions of are accompanied by large EP matrix 
elements. We have verified that the Q=(|,|,0)^ 
T\ (unstable) phonon (near K) causes large band 
shifts with atomic displacement (Sek/Su w 5 eV/A) 
near the FS necks, while for the stable T2 mode band 
shifts are no more than 5% of this value. Thus the fo- 



cusing of scattering processes is indeed coupled with 



TABLE I: From the calculated a 2 F(uj) at three pressures 
(GPa), the logarithmic and second moments of the fre- 
quency (K), the value of A, and T c (K) calculated using 
/i*=0.13. 



Pressure 






A 


T c 





209 


277 


0.40 


0.4 


10 


225 


301 


0.65 


5 


20 


81 


176 


3.1 


20 



large, polarization-dependent matrix elements. 

This focusing of EP coupling strength makes ac- 
curate evaluation of the total coupling strength A 
numerically taxing. The richness and strong in- 
dependence of the electron-phonon coupling that we 
have uncovered may explain the overestimates of T c 
in the previous work in Li, and may apply to the 
overestimates in boron [24|. It is clear however that it 
is EP coupling and not Coulomb interaction[25| that 
is responsible for the impressively high T c . Com- 
pressed Li thus has several similarities to MgB2 - 
very strong coupling to specific phonon modes, T c 
determined by a small fraction of phonons - but 
the physics is entirely different since there are no 
strong covalent bonds and it is low, not high, fre- 
quency modes that dominate the coupling. Com- 
pressed Li is yet another system that demonstrates 
that our understanding of superconductivity arising 
from "conventional" EP coupling is far from com- 
plete, with different systems continuing to unveil un- 
expectedly rich physics. 
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